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[02] GOVERNMENT RIGHTS 
5 [03] This invention is related to work performed in contract with the U.S. Government 
under DARPA ITO Contract #N66001-01-C-8042, "DARPA NeoCAD", and the 
U.S. Government may have certain rights in this invention. 

[04] PRIORITY CLAIM 
1 0 [05] This application claims the benefit of priority to provisional patent application 
60/418,044, filed October 12, 2002 with the U.S. Patent and Trademark Office, 
and titled "Method, Apparatus, and Computer Program Product for Simulation of 
Mixed-signal Circuits." 

15 [06] BACKGROUND OF THE INVENTION 

[07] (1) Technical Field 

[08] The present invention relates to techniques for system simulations. More 

specifically, the present invention relates to a technique for using wavelet operator 
to simulate the performance of mixed-signal circuits. 

20 

[09] (2) Discussion 

[10] Mixed-signal circuits, an example of which is the delta-sigma modulator, include 
both analog and digital functionality on the same chip and are difficult to simulate 
with conventional CAD software, such as SPICE or Simulink, for three main 
25 reasons: (1) they are described by a large number of equations; (2) the equations 

involve highly discontinuous non-linear operations at the clock period of the 
digital circuit; and (3) the equations are currently solved using slow, time- 
marching, algorithms (Runge-Kutta type). 
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[11] Recently, several approaches to fast simulation of mixed-signal circuits have been 
presented. Several are listed in the set of references below and are described here. 
Opal et al. presented a basic approach for circuits with a clock period in which 
linear differential equations are solved by one matrix multiply per clock cycle. In 
their method, a strong nonlinearity, such as the quantizer in a delta-sigma 
modulator, is simulated with a behavioral model at each clock period. Schreier 
and Zhang use a similar approach to construct recursion relations that update state 
variables of a delta-sigma modulator from time t to time t + T c> where T c 
represents the clock period. Cherry and Snelgrove compare three approaches: the 
recursion relation or direct integration approach, the time-marching method, and a 
z-domain extraction procedure, which were intended to combine the speed of the 
recurrence relations with the versatility of the time-marching method. Zhou et al., 
and Meliopoulos and Lee, have considered wavelet methods for use in general 
nonlinear circuit simulation and transient analysis. 

[12] There are several problems with the above-mentioned techniques. The matrix 
operator and direct integration methods are not sufficiently general to solve 
realistic circuit problems. The matrix in the matrix method is difficult to derive 
and is not small or sparse in general. The time-marching method is slow. The z- 
domain extraction procedure is difficult to generalize for real circuits. The 
wavelet techniques considered by Zhou et al., and Meliopoulos and Lee, could not 
be applied to mixed-signal circuits because of the number of wavelets required for 
time simulation over an entire time interval of 2 14 or more clock periods needed 
for analyzing mixed-signal circuits. 

[13] Accordingly, there exists a need in the art for a technique that overcomes the 
aforementioned limitations and that permits the use of wavelet operators to 
simulate the performance of systems in general, and mixed-signal circuits in 
particular. 
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[14] The following references are provided for additional, non-critical, information 

that may be of use to the reader (all are hereby incorporated herein by reference in 
their entirety): 

[15] [1] Y. Ahmed and A. Opal, "An efficient simulation method for oversampled 
delta-sigma modulators," Proceedings of the 37th Midwest Symposium on 
Circuits and Systems, vol 2 , 1994, pp. 1 164 -1 167. 

[16] [2] A. Opal, "Sampled data simulation of linear and nonlinear circuits," IEEE 

Trans. CAD Integrated Circuits and Syst., vol. 15, no. 3, Mar. 1996, pp. 295-307. 

[17] [3] Y. Dong and A. Opal, "Time-domain thermal noise simulation of switched 
capacitor circuits and delta-sigma modulators," IEEE Trans. CAD Integrated 
Circuits and Syst., vol. 19, no. 4, Apr. 2000, pp. 473-481. 

[18] [4] R. Schreier and B. Zhang, "Delta-Sigma modulators employing continuous- 
time circuitry," IEEE Trans. Circuits and Syst. I, vol. 43, no. 4, Apr. 1996, pp. 
324-332. 

[19] [5] J. A. Cherry and W. M. Snelgrove, "Approaches to simulating continuous- 
time delta-sigma modulators," Proceedings of the 1998 IEEE International 
Symposium on Circuits and Systems, ISCAS '98, vol. 1, 1998, pp. 587 -590. 

[20] [6] D. Zhou, W. Cai, and Wu Zhang, "An adaptive wavelet method for nonlinear 
circuit simulation," IEEE Trans. Circuits and Syst. I, vol. 46, no. 8, Aug. 1999, 
pp. 931-938. 

[21] [7] A. P. S. Meliopoulos and C.-H. Lee, "An alternative method for transient 

analysis via wavelets," IEEE Trans. Power Delivery, vol. 15, no. 1, Jan. 2000, pp. 
1 14-121 . A. W. Galli, "Discussion of 'An alternative method for transient 
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analysis via wavelets,' " Ibid, no. 4, Oct. 2000, p. 1326. A. P. S. Meliopoulos and 
C.-H. Lee, "Closure to discussion of 'An alternative method for transient analysis 
via wavelets,' " Ibid, no. 4, Oct. 2000, pp. 1326-1327. 



5 [22] [8] M. Unser and A. Aldroubi, Polynomial Splines and Wavelets — A Signal 

Processing Perspective, in Wavelets: A Tutorial in Theory and Applications, ed. 
by C. K. Chui, Academic Press, New York, 1992, pp. 91-122. 

[23] [9] H. L. Resnikoff and R. O. Wells, Jr., Wavelet Analysis, Springer, New York, 
10 1998, pp. 236-265, 281-340. 

[24] [10] G. Raghavan, J. F. Jensen, J. Laskowski, M. Kardos, M. G. Case, M. 

Sokolich, and S. Thomas III, "Architecture, design, and test of continuous-time 
tunable intermediate- frequency bandpass delta-sigma modulators," IEEE J. Solid- 
15 State Circuits, vol. 36, no. 1, Jan. 2001, pp. 5- 13. 

[25] [1 1] T. Kaplan, P. Petre and G. C. Valley, "State variable simulation of non- 
idealities in a continuous time delta-sigma modulator," submitted to IEEE J. CAD 
Integrated Circuits and Syst., May 2002. 

20 

[26] [12] W. H. Press, S. A Teukolosky, W. T. Vetterling, and B. P. Flannery, 

Numerical Recipes in Fortran, the Art of Scientific Computing, Second Edition, 
Cambridge University Press, 1992, pp. 340-386. 

25 [27] [13] O.C. Zienkiewicz and R. L. Taylor, Finite Element Method: Volume 1 The 
Basis, Butterworth and Heinemann, Oxford, 2000. 

[28] [14] G. Beylkin, "On the representation of operators in bases of compactly 
supported wavelets", SLAM Journal of Numerical analysis, vol. 6, no. 6, 
30 December 1992, pp. 1716-1740. 
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[29] SUMMARY OF THE INVENTION 
[30] The present invention comprises a method, an apparatus, and a computer program 
product for simulating a mixed-signal system. In one aspect, the invention 



15 [31] In another aspect, the system modeled is an electrical circuit, and the circuit may, 
more specifically, be a delta-sigma modulator. 

[32] In yet another aspect, in the generating operation, the matrix-based wavelet 
operator is developed by a wavelet-Galerkin method. 



[33] In a further aspect, in the generating operation, the matrix-based wavelet operator 
is developed directly from a system diagram or from equations that describe the 
system. 

25 [34] In a yet further aspect, in the operation of selecting, the number of wavelets is 
selected independently for each iteration of the acts of the method. 

[35] In a still further aspect, in the selecting operation, the set of wavelet basis 

functions is selected independently for each iteration of the acts of the method. 



10 



5 



comprises operations of generating a matrix-based wavelet operator 
representation of equations characterizing a system, with the matrix-based 
wavelet operator representation including wavelet connection coefficients. The 
next operation is selecting a number of wavelets and a set of wavelet basis 
functions with which to represent a performance of the system, whereby the 
wavelet operator, the number of wavelets and the set of wavelet basis functions 
represent a wavelet model of the system. The next operation is iteratively 
applying the wavelet model over a series of clock cycles to develop a behavioral 
model of the system. 



20 



30 
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In a still further aspect, the invention further comprises operations of receiving a 
specification for a system model and outputting the behavioral model of the 
system. 

One skilled in the art will appreciate that the operations described here have 
parallels in acts of a method, means of an apparatus, and means of a computer 
program product. In particular, the means of an apparatus and of a computer 
program product may be in the form of either software or a hard-coded instruction 
set. 

In addition, one skilled in the art will appreciate that a variety of additional 
aspects may be derived from combinations of the above aspects, and that these 
combinations are considered within the scope of the present invention. 

[39] BRIEF DESCRIPTION OF THE DRAWINGS 
The objects, features and advantages of the present invention will be apparent 
from the following detailed descriptions of the preferred aspect of the invention in 
conjunction with reference to the following drawings. 

FIG. 1 is a block diagram of a computer system that may be used in conjunction 
with the present invention; 

FIG. 2 is an illustrative diagram of a computer program product embodiment of 
the present invention; 

FIG. 3 is a flowchart depicting the general operations of the present invention; 

FIG. 4 is a circuit diagram of a continuous-time delta-sigma modulator used as an 
example of a system to which the present invention may be applied for analysis; 
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[45] FIG. 5(a) is a graph depicting the frequency response of the delta-sigma 

modulator to a 980 MHz cosine input signal, analyzed using 8, 16, and 32 Haar 
wavelets per state variable, and compared with a time marching solver; 

[46] FIG. 5(b) is an expanded view of a portion of the graph shown in FIG. 5(a); 

[47] FIG. 6(a) is a graph depicting the frequency response of the delta-sigma 

modulator to a 980 MHz cosine input signal, analyzed using a fast wavelet 
nonlinear solver, and compared with a time marching solver; and 

[48] FIG. 6(b) is an expanded view of a portion of the graph shown in FIG. 6(a). 

[49] DETAILED DESCRIPTION OF PREFERRED ASPECTS 
[50] The present invention relates to techniques for system simulations. More 
specifically, the present invention relates to a technique for using wavelet 
operators to simulate the performance of mixed-signal circuits. The following 
description, taken in conjunction with the referenced drawings, is presented to 
enable one of ordinary skill in the art to make and use the invention and to 
incorporate it in the context of particular applications. Various modifications, as 
well as a variety of uses in different applications, will be readily apparent to those 
skilled in the art, and the general principles defined herein, may be applied to a 
wide range of aspects. Thus, the present invention is not intended to be limited to 
the aspects presented, but is to be accorded the widest scope consistent with the 
principles and novel features disclosed herein. Furthermore, it should be noted 
that unless explicitly stated otherwise, the figures included herein are illustrated 
diagrammatically and without any specific scale, as they are provided as 
qualitative illustrations of the concept of the present invention. 

[51] In order to provide a working frame of reference, first a glossary of terms used in 
the description and claims is given as a central resource for the reader. Next, a 
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discussion of various physical aspects of the present invention is provided. 
Finally, a discussion is provided to give an understanding of the specific details. 

(1) Glossary 

Before describing the specific details of the present invention, a centralized 
location is provided in which various terms used herein and in the claims are 
defined. The glossary provided is intended to provide the reader with a general 
understanding of the intended meaning of the terms, but is not intended to convey 
the entire scope of each term. Rather, the glossary is intended to supplement the 
rest of the specification in more accurately explaining the terms used. 

Means - The term "means" as used with respect to this invention generally 
indicates a set of operations to be performed on a computer, and may represent 
pieces of a whole program or individual, separable, software modules. Non- 
limiting examples of "means" include computer program code (source or object 
code) and "hard-coded" electronics (i.e. computer operations coded into a 
computer chip). The "means" may be stored in the memory of a computer or on a 
computer readable medium. 

(2) Physical Aspects 

The present invention has three principal "physical" aspects. The first is a 
system, typically in the form of a computer system operating software or in the 
form of a "hard-coded" instruction set. This system may either serve as the basis 
for a stand-alone system analysis package or may be integrated with other 
computer-aided design (CAD) and simulation programs. The second physical 
aspect is a method, typically in the form of software, operated using a data 
processing system (computer). The third principal physical aspect is a computer 
program product. The computer program product generally represents computer 
readable code stored on a computer readable medium such as an optical storage 
device, e.g., a compact disc (CD) or digital versatile disc (DVD), or a magnetic 
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storage device such as a floppy disk or magnetic tape. Other, non-limiting 
examples of computer readable media include hard disks, read only memory 
(ROM), and flash-type memories. These aspects will be described in more detail 
below. 

A block diagram depicting the components of a computer system used in the 
present invention is provided in FIG. 1. The data processing system 100 
comprises an input 102 for receiving specifications of a system - this may be in 
the form of a system diagram or in the form of equations representing a system. 
Note that the input 102 may include multiple "ports". Typically, input is received 
through standard input/output devices, non-limiting examples of which include a 
keyboard, mouse, hard drive, floppy drive, network connection, etc. The output 
104 is connected with the processor for providing information regarding the 
analysis for display to a user via devices such as printers and computer monitors. 
Output may also be provided to other devices or other programs, e.g. to other 
software modules, for use therein. The input 102 and the output 104 are both 
coupled with a processor 106, which may be a general-purpose computer 
processor or a specialized processor designed specifically for use with the present 
invention. The processor 106 is coupled with a memory 108 to permit storage of 
data and software to be manipulated by commands to the processor. 

An illustrative diagram of a computer program product embodying the present 
invention is depicted in FIG. 2. The computer program product 200 is depicted as 
an optical disk such as a CD or DVD. However, as mentioned previously, the 
computer program product generally represents computer readable code stored on 
any compatible computer readable medium. 

(3) Introduction 

The present invention provides a technique that can be used to simulate the 
performance of a wide variety of systems, non- limiting examples of which 
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include mechanical and electrical systems. Although it is readily applicable to the 
analysis of a broad range of system types (mechanical, electrical, etc.), the present 
invention will be discussed herein in the context of calculating the temporal 
response of mixed-signal circuits. In particular, such a system has application to 
the next generation of computer aided design (CAD) software for circuit design 
and analysis. It is contemplated that one of skill in the art would readily be able 
to adapt the techniques presented herein to other systems without difficulty. 



[61] The present invention provides a technique in which a system simulation is 
10 broken up into clock periods, and the calculation within each clock period, which 

is weakly non-linear, is performed by matrix multiplication. Also, the technique 
allows the discontinuous non-linearity to be applied separately at the end of the 
clock period. The entire simulation can be performed by repeating this process 
for the required number of clock periods. Since the technique presented is much 
1 5 faster than time-marching algorithms that require small time steps to simulate the 

discontinuous non-linearity, a much larger number of equations can be simulated 
in reasonable computer run times. 



[62] In one aspect, the present invention uses the wavelet-Galerkin method to derive a 
20 matrix representative of the system to be analyzed. This matrix, which is used 

during each clock period of the analysis, is derived from a set of differential 
equations that describe the system. This matrix can also be directly derived from 
a circuit-type schematic of the system. Techniques are also described for 
optimizing the representation of the system. Examples are provided to clearly 
25 explain how the invention may be applied. Specifically, for the purpose of 

explanation, a 4 th order delta-sigma modulator is used as an example, and the 
results are in excellent agreement with conventional calculations. The invention 
may be applied to any time-based mixed-signal system that operates on a clock- 
basis with strong or discontinuous non-linearity at the clock period and that can 
30 be described by conventional linear or non-linear, integral or differential, 
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equations for the state variables (e.g., current, voltage, or temperature) during the 
clock period. 

The invention provides for faster simulations of mixed-signal systems than the 
simulations obtained with conventional techniques. Furthermore, the speed 
advantage permits more calculations to be performed during a given time period 
for system optimization. Thus, this invention can be used to solve a more 
complex set of equations in the same amount of time so that low-level (transistor- 
level, in the case of electrical circuits) simulations may be performed for systems 
that were previously only capable of being simulated at the behavioral level. The 
use of wavelet expansion techniques to simulate system performance permits 
higher resolution simulations for the same run time, because the wavelets can be 
chosen to fit the details of the circuit. 

Each clock period in a simulation requires one wavelet transform of the input and 
any applicable feedback signals in the system, but the transform can be performed 
using a fast wavelet transform technique. The wavelet-Galerkin matrix operator, 
which is used during a clock period, can be derived directly from the circuit 
layout. This feature makes the technique of the present invention attractive for 
applications such as CAD. 

Non-linearities within the clock period can be simulated with the techniques 
presented by using any one of many well-known iterative solving techniques with 
the fast wavelet transform and the collocation point method. Higher-order partial 
differential equations for state variables can be simulated with this technique by 
setting the first derivative of the state variable equal to a new state variable and 
adding an equation for this new state variable; or directly in the wavelet-Galerkin 
approach by using a wavelet basis for which the connection coefficients for the 
higher-order derivatives exists. 
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The general operations performed by the present invention are depicted in the 
flowchart shown in FIG. 3. In the first operation, a matrix-based wavelet operator 
representation of equations characterizing a system is generated 300, with a 
matrix including wavelet connection coefficients. Next, a number of wavelets 
and a set of wavelet basis functions are selected 302 for representing the 
performance of the system. The wavelet operator, the number of wavelets, and 
the set of wavelet basis functions represent a wavelet model of the system. This 
wavelet model is applied 304 over a series of clock cycles to develop a behavioral 
model of the circuit (the repetitive application is represented by the dashed line in 
the diagram). 

The approach presented herein will now be explained in greater detail by use of 
an example. The example presented is an electrical circuit. It is important to 
note, however, that the present invention is not limited either to this specific 
electrical circuit, or to electrical circuits in general, and that it may be applied to a 
broad range of systems. 

(4) Exemplary Description 

a. Delta-Sigma Modulator Specification 

The present invention will now be described in detail with reference to a specific 
example of a mixed-signal system, specifically, a 4 th order continuous-time delta- 
sigma modulator (a mixed-signal electrical circuit). For more background 
regarding the design and implementation of continuous-time delta-sigma 
modulators, the reader is referred to the reference by G. Raghavan et al., cited in 
the set of references above. Except as provided herein, specific design issues 
related to delta-sigma modulators are not considered necessary to understand or 
practice the present invention. 
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[71] A circuit diagram of a continuous-time delta-sigma modulator used in this 
example is presented in FIG. 4. As shown, the delta-sigma modulator 400 
comprises an input 402 leading to the first of a set of transconductance amplifiers 
404 (go-s) that are arranged as shown. The delta-sigma modulator 400 also 
includes a plurality of adders 406, transimpedance integrators 408 (1/sCj), a single 
quantizer 410, a set of delays 412, a set of digital-to-analog converters (DACs) 
414 (#6-9), and a voltage amplifier 416 (0fb). An output 418 of the quantizer 410 
is simulated in this example. Note that the symbols in the figure are summarized 
as follows: 

[72] g\ represent gain cells; 

[73] + represent adders; 

[74] Q represent integrators; 

[75] At represent time delays; 

[76] Q represents the quantizer; and 

[77] u(t) represents the input to the delta-sigma modulator. 

[78] The equations for the delta-sigma modulator of FIG. 4 are as follows: 



[79] 



at 



= go u{t) - g 2 x 2 (t) - g 6 y(t - At) 



(1) 



[80] 



dx 2 (t) 
dt 



= gixi(t)-gi y(t-kt) 



(2) 
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c$ — — = g$ x 2 (t) - g s x 4 (t) - gs yit - A/) 

[81] dt ; (3) 



dx 4 (t) 

[82] dt ; and (4) 

5 

[83] y = sgn[jc 4 (f„) + Q^y{t n A)\\ where (5) 

[84] jtj are the voltages following the integrators Q; 

10 [85] yis the output of the quantizer; and 

[86] the other symbols are as mentioned previously. Also, the clock period is T c and 
t„= n T c is the time at the end of the « th clock cycle. 

15 [87] b. Derivation of the Wavelet Operator 

[88] Given the delta-sigma modulator definition just provided, the next step is to 
convert the equations (l)-(5) into a matrix equation. To do so, the wavelet- 
Galerkin method presented in H. L. Resnikoff and R. O. Wells, Jr., Wavelet 
Analysis, Springer, New York, 1998, pp. 236-265, 281-340 (incorporated by 

20 reference herein in its entirety) is used. Note that the use of the wavelet-Galerkin 

method is for exemplary purposes only, and is not intended to limit the scope of 
the present invention. 

[89] The bounded domain in this case is a single clock period, and the wavelet- 
25 Galerkin method is repetitively applied for the number of clock periods in the full 

simulation. Since the output 418 y(f n ) is held fixed for t n < t < t n + T C9 both the 
output 418 y(t) and the input 402 u(t) are known functions during the clock period 
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and remain on the right hand side of equations (l)-(4). These functions plus the 
boundary conditions for the state variables are expanded in a set of wavelet basis 
functions. Then, the state variables are expanded in the same set of wavelet basis 
functions that enable the wavelet transform of equations (l)-(4) to be put in a 
matrix representation: 



[90] 
[91] 



10 





8ll 


0 


0 


-gil 


c 2 r 


0 


0 


0 




CsT 


Ssl 


0 


0 


-84 1 


c 4 r 



X = R 



, where 



(6) 



T represents a matrix of the wavelet connection coefficients that relates the 
wavelet expansion of the derivative of a function to the expansion of the function 
itself (equivalently, the connection coefficients relate the derivative of the scaling 
and wavelet functions to the functions themselves); 



15 



[92] / is the identity matrix; and 

[93] X is the voltage vector corresponding to the wavelet expansion coefficients of X\ 
through X4 and R to the expansion coefficients of the right hand sides of equations 
(l)-(4). 



20 [94] 



25 



If m wavelets basis functions are used per state variable, then the dimensions of 
T 5 J , and 0 are m x m, while X and R are 4 m long. It is evident that the matrix in 
Equation (6) is quite sparse (r and / are also sparse sub-matrices), and the other 0 
matrices on the left of Equation (6) reflect the fact that there are relatively few 
connections between the circuit nodes shown in FIG. 4. It should be noted that 
Equation (6) does not depend on the number of wavelet basis functions or the type 
of wavelet basis functions. The invention requires solely that the wavelet basis 
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functions chosen for analysis have a connection coefficient matrix. Note that the 
article by H. L. Resnikoff et al, just mentioned, provides explicit formulas for 
connection coefficients for wavelets of genus 1 to 6 and suggests, although it is 
not mathematically proven, that the connection coefficients may also exist for 
other wavelets. It is also noteworthy that Equation (6) does not solve equations 
(l)-(5) simultaneously, but only equations (l)-(4) within a clock period during 
which y is a known function from previous clock cycles. The technique presented 
here includes solving Equation (6) by inversion of the matrix on the left-hand side 
and multiplication of the vector R by this matrix. This yields the wavelet 
expansion coefficients for the state variables. One final inverse wavelet transform 
could be used to obtain the full time-dependent state variables within the clock 
period; but since only the values of the state variables are required at the clock 
period, a full wavelet transform is not necessary for this example. After the 
wavelet expansion coefficients for the state variables have been obtained, the 
quantizer 410 (see FIG. 4 and Equation (5)), input signal 402, and boundary 
conditions are updated, and the process is repeated for the next clock period. A 
complete simulation of the delta- sigma modulator typically requires 2 14 or more 
clock periods that corresponds to 2 14 matrix multiplications plus 2 14 wavelet 
transforms. For each choice of the wavelet basis and the number of wavelets, the 
matrix on the left hand side of Equation (6) must be inverted once, but this is 
unimportant because the time is amortized over 2 14 or more clock cycles. For 
computational efficiency, the wavelet transform may be performed with the fast 
wavelet transform algorithm. Further, the number of wavelet basis functions need 
not be constant from one clock period to another. Larger numbers of basis 
functions might be used when the input signal varies more rapidly for instance. 
Similarly, the type of wavelet basis function need not be the same from clock 
period to clock period. For example, one might use a Haar basis for one clock 
period and a Daubechies basis for another along with an intelligent adaptive 
algorithm to choose from a library of pre-computed matrices. 
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c. Derivation of the Wavelet-Galerkin Matrix from the Circuit Diagram 
Note that the wavelet-Galerkin matrix on the left hand side of Equation (6) may 
be written down by inspection from the circuit diagram in FIG. 4. The procedure 
used is as follows. Start at the first integrator C\. Since the output of this 
integrator is jci, it gives rise to a term C\T in the first row of the matrix. The input 
to the integrator consists of go u(t\ g 6 y(t n -At) and g2X 2 . The wavelet transform of 
the sum of the first two terms is the first element of the vector R. The term #2*2 
gives the sub-matrix gj in the matrix. Similar analysis at the other integrators 
provides the other rows of the matrix. 

d. Reduction to Practice 

The present invention has been reduced to practice with 8, 16, and 32 Haar 
wavelets per state variable Xi and these results have been compared to those from 
a time marching code, which has been independently benchmarked to Simulink 
calculations (Simulink is a product of The Math Works Inc., 3 Apple Hill Drive, 
Natick, MA 01760-2098, United States). The result of these comparisons are 
shown in FIG. 5(a) and FIG. 5(b). FIG. 5(a) shows the absolute value of the 
discrete Fourier transform of the output of the delta-sigma modulator multiplied 
by the Kaiser Window as a function of frequency for a 4 th order delta-sigma 
modulator signal over the full range of frequency from 0 to 2 GHz, and FIG. 5(b) 
shows the same with an expanded frequency scale showing a sharp notch in 
modulator frequency response at 1.016 GHz. These calculations use the 
parameters of the delta-sigma modulator shown in FIG. 4, and assume a 980 MHz 
cosine signal as input, a clock period of 250 ns, and a DAC delay of TJ2. The 
output 418 from the quantizer 4 10, y, consists of atableof l's and-l's. The 
spectrum is obtained by multiplication by a Kaiser window function and taking 
the discrete Fourier transform to obtain a spectral representation of the output 
signal from the circuit. Note that on the full-scale plot in FIG. 5(a), the results 
(signal at 980 MHz and average noise level) appear nearly identical regardless of 
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the number of wavelets used per clock cycle, with minor differences in the noise 
spikes. In the blown-up version shown in FIG. 5(b), it is shown that the wavelet 
operator code with 8 and 16 Haar wavelets does not predict the sharp "notch" to 
the right of the signal. This notch is a very sensitive indicator of the quality of the 
simulation. 

[99] e. Extension to Circuit with Nonlinear Elements 

[100] The fast wavelet operator method may be extended to circuits that contain 
nonlinear elements. To illustrate use of the wavelet operator method with a 
nonlinearity in the delta-sigma modulator, it is assumed that gain cell g\ 404 (see 
FIG. 4) saturates according to the semi-empirical expression (the algorithm is 
independent of the detailed form of the nonlinearity), 



[102] where v is the voltage (jci for cell g\) and v s is a constant that characterizes the 
saturation. The right hand side of the Substitution (7) has a typical bell shaped 
dependence, e.g., it equals 1 for v « v s and 0 for v » v s . Substitution (7) is 
obtained by fitting an I-V curve for a typical gain cell. 

[103] Substitution (7) can be used directly in Equation (2) in the time marching code 
without further modification. These results, shown in FIGs. 5(a) and 5(b), 
provide the baseline against which to compare the wavelet results. Specifically, 
FIG. 5(a) and 5(b) depict noise power/bin as a function of frequency for the 
exemplary 4 th order delta-sigma modulator. The wavelet code for this non- 
linearity makes use of an iterative wavelet solver. Examples of iterative solvers 
are commonly known to those of skill in the art, and, for the benefit of the reader, 
several examples are provided in the references in the background. To use the 
one method, the Equation (2) is rewritten into the following form: 



[101] 



g-> g 



( 1+ i) sech2 [^( 1+ A)], 



(7) 
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dx t (t) 
[104] dt 
. (8) 



-*7^-A#)- ft x|(0 + ft [l + Jsech 2 [— — j^l + —j-jj^d) 



5 [105] which reduces to Equation (2) forx\(t) « v s . As in the linear case, Haar wavelets 
are used to illustrate the method, but the method can be used with arbitrary 
wavelets. The first iteration at each clock cycle is obtained as described above 
with x\(t) set to 0 on the right hand side of Equation (8). The second iteration is 
obtained by using the same wavelet operator but with an added source term on the 

10 right hand side corresponding to the saturation term evaluated with the value of 

x\(t) obtained in the first iteration. Since the only knowledge of x\(t) is its 
wavelet coefficients, this might seem to be a computationally intensive task. 
Fortunately, existing wavelet computation technology is available for this 
purpose. A combination of collocation points and fast wavelet coefficient 

15 computations is used for this operation. Essentially, for N wavelets, the nonlinear 

evaluation needs to be performed only N times at the collocation points, and 
calculation of the new coefficients of the right hand side is completed by the fast 
wavelet transform. 



20 [106] Since the delta- si gma modulator functions poorly, if at all, with large nonlinear 
gain saturation, for practical cases of interest, the number of iterations required is 
less than 10 and often less than 5 for good accuracy. It is obviously advantageous 
to adaptively control the number of iterations in each clock cycle since the value 
of x\(t) is rarely large enough to require much iteration. For mixed signal circuits 

25 in which the non-linearity is large within the clock period, it may be necessary to 

use other iterative solvers. All classes of iterative solves can be used very 
efficiently because the linear equations are solved with one matrix multiply and 
the right hand sides are evaluated with the fast wavelet transform. 
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[107] The wavelet operator and the time marching results for the 4 th order delta-sigma 
modulator with gain cell g\ 404 saturated are shown in FIGs. 6(a) and 6(b). These 
results are obtained with 8 Haar wavelets and correspond to the case shown in 
FIGs. 5(a) and 5(b), but with saturation voltages of v s = 1.72 and 2.35 V (for 
FIGs. 6(a) and 6(b), respectively. 

[108] f. Extensions of the Wavelet-Galerkin Method 

[109] Several extensions of the wavelet-Galerkin algorithm will be readily ascertainable 
by one skilled in the art of wavelet analysis. First, the wavelet basis can be 
chosen to match properties of the circuit. This will help to reduce the size of the 
matrix. Second, the same wavelet basis need not be used for each clock cycle. 
An adaptive procedure can be devised so that the wavelet basis changes 
depending on properties of the input signal, the feedback DAC or other circuit 
non-idealities. Note that the calculations presented herein were illustrated with a 
very simple model of the DAC, but the method is not limited to a specific 
description of the DAC. The properties of the DAC may, however, influence the 
choice of a wavelet basis. Finally, the number of wavelets used per clock period 
can vary adaptively depending on the input conditions, the DAC, and any other 
circuit non-idealities that may be present. It is believed that other extensions of 
the wavelet-Galerkin method disclosed in this work will be clear to those skilled 
in conventional Galerkin methods. 

[110] g. Extension to Other State Variables 

[111] Complete description of the behavior of a system (mixed signal circuit in the 
example herein) may require other state variables besides current and voltage. 
For example, in the delta-sigma modulator temperature may affect DAC response. 
It should be clear that the x\ in equations having the same form as Equations (1)- 
(5) can correspond to arbitrary state variables that can be described by first order 
differential equations. 
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[112] h. Extension to Higher-Order Differential and Integral Equations 

[113] The method is not limited to state variables that can be described by first order 
differential equations. Differential equations with n th order derivatives may be 
reduced to n first order equations by repeated application of the following 
technique illustrated for a second order equation. If the original second order 
differential equation is dfx/d^fldxfdtjcj], and it may be replaced two 
equivalent first order equations: 

[114] u = dx/dt and 

[115] d*/d/ 

[116] An alternative method for state variables described by higher order differential or 
integral equations is to use the wavelet-Galerkin method with a wavelet basis for 
which the connection coefficient matrix for higher order integration or 
differentiation exists. 
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